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We describe a formulation to deduce the phase shifts, which determine the ground-state 
properties of interacting quantum-dot systems with the inversion symmetry, from the fixed- 
point eigenvalues of the numerical renormalization group (NRG) . Our approach does not assume 
the specific form of the Hamiltonian nor the electron-hole symmetry, and it is applicable to a 
wide class of quantum impurities connected to noninteracting leads. We apply the method to 
a triple dot which is described by a three-site Hubbard chain connected to two noninteracting 
leads, and calculate the dc conductance away from half-filling. The conductance shows the 
typical Kondo plateaus of the Unitary limit g ~ 2e 2 /h in some regions of the onsite energy 
td, at which the total number of electrons 7V e i in the three dots is odd, i.e., N e \ ~ 1, 3 and 
5. In contrast, the conductance shows a wide minimum in the regions of ed corresponding to 
even number of electrons N e i ~ 2 and 4. We also discuss the parallel conductance of the triple 
dot connected transversely to four leads, and show that it can be deduced from the two phase 
shifts defined in the two-lead case. 
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1. Introduction 

The Kondo effect in quantum dots is a subject of much 
current interest, 1,2 and experimental developments in 
this decade make it possible to examine the interplay of 
various effects which have previously only been studied 
in different fields of physics. 3,4 For instance, in quantum 
dots, the interplay of the Aharanov-Bohm, Fano, Joseph- 
son, and Kondo effects under equilibrium and nonequi- 
librium situations have been studied intensively. 5-7 

To study the low-temperatures properties of the sys- 
tems showing the Kondo behavior, reliable theoretical 
approaches are required. The Wilson numerical renor- 
malization group (NRG) method 8-10 has been used suc- 
cessfully for single and double quantum dots. 11,12 In a 
previous work, we have applied the NRG method to a fi- 
nite Hubbard chain of finite size Nc connected to nonin- 
teracting leads, and have studied the ground-state prop- 
erties at half- filling. 13 The results obtained for Nc = 3 
and 4 show that the low-lying eigenstates have one-to- 
one correspondence with the free quasi-particle excita- 
tions of a local Fermi liquid. It enables us to determine 
the transport coefficients from the fixed-point Hamilto- 
nian. Furthermore, it enables us to deduce the character- 
istic parameters such as Tk and Wilson ratio R from the 
fixed-point eigenvalues, and the calculations have been 
carried out for the Anderson impurity model. 14, 15 The 
purpose of this paper is to provide an extended formula- 
tion to deduce the conductance away from half-filling. 

The ground-state properties of quantum dots con- 
nected to two noninteracting leads are determined by the 
two phase shifts if the systems have an inversion sym- 
metry. Kawabata described the outline of this feature 
assuming that the low-energy properties are determined 
by the quasi-particles of a local Fermi liquid, 16 and dis- 
cussed qualitatively the role of the Kondo resonance on 
the tunneling through a quantum dot. In this paper we 



describe the formulation to deduce these two phase shifts 
from the fixed-point eigenvalues of NRG. The formula- 
tion along these lines is well known for some particu- 
lar cases, and has been applied to the single and double 
quantum dots. 5, 11,12 However, to our knowledge, a gen- 
eral description has not been given explicitly so far. Our 
formulation does not assume a specific form of the Hamil- 
tonian, and thus it is applicable to a wide class of the 
models for quantum impurities. We apply the method to 
a three-site Hubbard chain. It is a model for a triple dot, 
and the low-energy properties of this system away from 
half-filling have not been clarified with reliable methods 
yet. We calculate the dc conductance g away from half- 
filling as a function of the onsite energy ta which can 
be controlled by the gate voltage. The results show the 
typical Kondo plateaus of the Unitary limit g ~ 2e 2 /h 
when the average number of electrons N c \ in the chain is 
odd, N c \ ~ 1, 3, and 5. In contrast, for some finite ranges 
of ed corresponding to even number of electrons N c i ~ 2 
and 4, the conductance shows wide minima. We also dis- 
cuss the parallel conductance of the triple dot connected 
transversely to four leads. 

In §2, we introduce a model for interacting electrons 
connected to reservoirs in a general form, and describe 
the ground-state properties in terms of the phase shifts. 
In §3, we describe the formulation to deduce the phase 
shifts from the fixed-point Hamiltonian of NRG. In §4, 
we apply the formulation to the transport through the 
triple dot. In §5, discussion and summary are given. 

2. Model and Formulation 

We start with a system consisting of a finite central 
region (C) and two leads, on the left (L) and the right 
(R), as illustrated in Fig. 1. The central region consists 
of Nc quantized levels, and the interaction Uj i j 3 -j 2 j 1 is 
switched on for the electrons staying in this region. Each 
of the two leads has a continuous energy spectrum. The 



2 J. Phys. Soc. Jpn. 



Full Paper 



A. Oguri, Y. Nisikawa, and A. C. Hewson 



Left lead 



Sample 



Right lead 




Fig. 1. Schematic picture of a series connection 



complete Hamiltonian is given by 
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where di a annihilates an electron with spin a at site i. 
In the lead at v (= L, R), the operator c\ U(y creates 
an electron with energy tku corresponding to an one- 
particle state 4>ki/(r). The tunneling matrix elements v L 
and v R connect the central region and two leads. At the 
interfaces, a linear combination of the conduction elec- 
trons t/Vo- = J2k c kua ^kvirv) mixed with the electrons 
at i = 1 or Nq, where r„ denotes the position at the in- 
terface in the lead side. We assume that the hopping 
matrix elements ty and v u are real, and the interac- 
tion has the time-reversal symmetry: U^i\ is real and 

^43;21 = C^34;12 = f/l2;34 = ^42;31 = ^13;24- We will be 

using units h — 1 unless otherwise noted. 

We use the Green's function for the interacting region 
(ij G C) defined by 



Gij (kj r , 



dr (T T d ia (r)dl(0) 



(6) 



where f3 = 1/T, d ja {r) 



djfjG 



-tH 



and (O) 



Tr [e @ n O] /Tre ^ . The corresponding retarded and 
advanced functions are given by (^.-(w) = Gij (to ±iS) 
via the analytic continuation. The Dyson equation for 
the Green's function G(z) = {Gij(z)} can be rewritten 
in an Nc x Nc matrix form 



{G(z)}- 1 = zl 
where H° c = {-%}, 

v 2 L g L (z) 
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and = {Sy(z)} is the self-energy due to the inter- 
electron interaction Tic- In the matrix Vmw(z), the 
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Fig. 2. Schematic picture of a parallel connection. 



two non-zero elements are determined by the Green's 
function at interface of the isolated lead g^(w) = 

J2k |^ CT ( r ^)| /( w " e fc^ + * 0+ )- We assume, for simplic- 
ity, that the density of states p v is a constant for small u. 
Then, the energy scale of the level-broadening becomes 
T,y = np u vl, and g+(w) = -i7rp„. 

2. 1 Ground state properties 

If the self-energy shows a property ImS + (0) = at 
T = 0, then the damping of the single-particle excitations 
at the Fermi level vanishes. In this case, the renormal- 
ization hopping matrix elements H 'q = { 

by 

off 



Uj } defined 



H 



c 



H 



c 



ReS+(0) 



(9) 



play a central role on the ground-state properties. The 
value of the Green's function at the Fermi level is given 
by {G + (0)}^ = K(0) - V+ x (0) with 



LOl H 



efi 



and the scattering coefficients at T 
by the determinant 



(10) 

are determined 



det{G + (0)} 1 = \-T L T R detK^ cNc {0)+detK(0) 



i T L detail (0) 



Tr dct K n c n c (0) 

(11) 



Here Kyfi) is an (N c - 1) X (N c - 1) derived from 
K(0) by deleting the i-th row and j-th column. Simi- 
larly, K^ cNc {0) is an {N c - 2) x {N c - 2) matrix ob- 
tained from K(0) by deleting the first and Nc-th rows, 
and the first and iVc-th columns. The dc conductance at 
T = can be expressed, using an inter-boundary Green's 



.9 S . 



2e^ 



^4r fl r L |G+ cl (o)| 



(12) 



Similarly, at T — the charge displacement iV e i caused 
by the hybridizations and Tr can be expressed, using 
the Friedel sum rule, 20 ' 21 as 



N cl = — log [det S] , 

7T1 



(13) 
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(14) 



In the case of the constant density of states we are 
assuming, the charge displacement coincides with the 
total number of electrons in the central region, i.e., 
N c i = ££°i where n d .„ = d\ a d la . Note 

that the wavefunction renormalization due to the matrix 
1 — dTi/{duj) does not contribute to the conductance and 
charge displacement at T = 0, although it has an infor- 
mation about Tk and can be related to the asymptotic 
behavior of the low-lying excited states near the fixed- 
point of NRG. 14 ' 15 ' 22 

Using the above Green's functions obtained for the 
series connection, one can calculate the conductance for 
the parallel connection, for which the tunneling matrix 
elements are given by vl/V% and for the leads 

labeled by "1" and "2", respectively as shown in Fig. 
2. Owing to this special symmetry between the left and 
right leads, the interacting site at i = 1 (i = Nc) is 
coupled only to an even combination of the states from 
the left and right leads with the label "1" ("2"). Thus, 
the Green's functions in the central region G(z) coincide 
with that for the series connection, so that a ground-state 
average such as iV e i becomes the same with that for the 
series connection. The parallel conductance, however, is 
different because two conducting channels contribute to 
the total current flowing in the horizontal direction, and 
it can be expressed in the form 



= — 

■^parallel h 



GtM\ +T% 



+ 2T L T R \G+(0) 



(15) 



This expression is for the current flowing from left to 
right (or right to left), and has been derived from the 
Kubo formula using a multi-channel formulation for the 
dc conductance. 21 

2.2 Inversion symmetric case 

We now consider the case the system has an inversion 
symmetry; T L = T R (= T), g L = g R (= g), and t itj = 
tN c -i+\,Nc-j+i- Due to this symmetry, the eigenstates 
of the Hamiltonian H can be classified according to the 
parity, and the matrix Dyson equation given in eq. (7) 
can be separated into the two subspaces corresponding 
to the even and odd orbitals; 



b 3,<? 
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where j — 1, 2, . . . , Nc/2 for even Nq, and j — 
1,2,..., (N c -l)/2 for odd N c . Note that there is one ex- 
tra unpaired orbitals a(jv c +i)/2,cr = ^(jv c +i)/2,ct f° r odd 
Nc- Thus the Green's functions for these orbitals a^ a 
and bj !<7 can be written separately in the forms 

{G+^O)}" 1 = K°™ (0) - V e ™(0) , (18) 



{G+ dd (0)} 1 = 2T> dd (0)-K? d (0), 



(19) 



at T = and u = 0. The matrices K-° von (0) and K odd (0) 
are derived from K(0) in eq. (10) via the transformation 
cqs. (16) and (17). Since only the even and odd orbitals 
with the label j = 1, i.e., ai ;Cr and &i ;Cr , have a finite 
tunneling- matrix element to the lead with the same par- 
ity, the mixing term defined in eq. (8) is transformed as 



V 7 . M = 












(20) 



even 
mix 



for 7 = "even" and "odd" . Note that the size of V 
and that of are identical for even Nc, whereas the 

size of these two matrices are different for odd Nc be- 
cause the one unpaired orbitals exists. The determinant 
given in eq. (11) can be factorized, as 

det{G+(0)} _1 = det{G+ 



1 (0)}- 1 det{Gf+ d (0)}- ] 
[|[dctK 7 (0) + iT dctK^O) 



(21) 



where the product runs over the two partial waves 7 = 
"even" and "odd". The matrices K^ cn (0) and K° dd (0) 
are derived respectively from i<: ovon (0) and K odd (0) by 
deleting the first row and column corresponding to the 
orbitals ai j(T and &i j(T , which are connected directly to 
nonintcracting leads. The scattering coefficients are de- 
termined by the ratios, n cven and K dd, defined by 



detK 7 (0) 

Ki-y 



(22) 



rdet-K^O) ' 

The value of K OV on and K 0( jd can be deduced from the 
fixed-point eigenvalues of NRG, as described in eq. (38) 
in §3. 

The phase shifts for the even and odd partial waves, 
(5 0VC n and <5 dd, can be defined using the retarded Green's 
functions for cii a and b l iy at the Fermi level, which can 
be expressed, using eq. (21), as 
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The Green's functions at the interfaces G^ c 1 (0) and 
G\ x (0) are given by the linear combinations 



G + Nc A0) = 
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Thus, the conductance defined in eq. (12) can be written 
in the form 



.9s 



sin 



(27) 



Similarly, the Friedel sum rule eq. (13) can be rewritten 
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Fig. 3. Schematic pictures of discretized Hamiltonian of NRG. 



in the form 

EEK-> = - (^ven + S odd j . (28) 
i=l a 

Furthermore, the parallel conductance for the current 
flowing in the horizontal direction in the geometry shown 
in Fig. 2 can be expressed using these two phase shifts, 
i.e., the inversion symmetry simplifies eq. (15) as 

2e 2 / \ 

^parallel = ~ ( sin2 ^even + sin 2 S odd J . (29) 

Thus the parallel conductance is a simple sum of the 
contributions of the even and odd channels. 11, 12 Alter- 
natively, it can be expressed in the form, <7 paral i e i = 
(2 e 2 /fe)2r{-ImG+(0)}. 

3. Fixed-point Hamiltonian and phase shifts 

In the NRG approach the conduction band is trans- 
formed into a linear chain as shown in Fig. 3 after carry- 
ing out a standard procedure of logarithmic discretiza- 
tion. 8-10 Then, a sequence of the Hamiltonian Hn is in- 
troduced, as 

H N = A( N -^(H%+H u c + H mix + Hl N J d ), (30) 

#mix = V J2 (/(U<A* + d l*f ,L*) 

(7 

+ « E ( /<U<W + 4 c , CT / ,^ ) - ( 31 ) 



H (N) - D 

n lead — u 



1 + 1/A 



N-l 

E EE^" /2 

i>=L,R a n=0 



X ( fn+l,i/a fn,va fn,va fn+l,va ) ' (^2) 

where D is the half-width of the conduction band. The 
hopping matrix elements v and are defined by 



2DYA K 



1 1 + 1/A 

2 1 - 1/A 



1/A 



n+l 



v/l - 1/A 2 "+Vl - 1/A 2 "+ 3 



(34) 



The factor A\ is required for comparing precisely the 
discretized Hamiltonian Hn with the original continuous 
model defined in eq. (I). 9 ' 23 In the discretized Hamilto- 
nian Hn, the matrix elements t and v are multiplied by 
A^- 1 )/ 2 ^ so that the original Hamiltonian Ti are recov- 
ered from A~( n ~ 1 ^ 2 Hn in the limit of N — > oo and 
A-» 1. 



In a wide class of impurity models described by eq. 
(1), the low- lying energy eigenvalues of Hn converge for 
large N to the spectrum which has one-to-one correspon- 
dence to the quasi-particles of a local Fermi liquid. 9 ' 10 ' 13 
It enables us to deduce the matrix elements of H b q from 
the NRG spectrum. At the fixed point, the low-energy 
spectrum of the many-body Hamiltonian Hn can be re- 
produced by the one-particle Hamiltonian consisting of 
H c g and the two finite leads; 



H W =A (N-l )/2 ^ HhS + Hmix + H] 



(AO 
lead 



(35) 



Here Hf 



7 , and tij is the renormal- 



ized matrix element defined in eq. (9). The Hamiltonian 
Hqp^ describes the free quasi-particles in the system con- 
sisting of Nc + 2(N + 1) sites, and the corresponding 
Green's function for the interacting sites can be written 
as an Nc x Nc matrix, 



{GqpM}- 1 EE 



A (JV-l)/2 



wA" 



-{N-l)/2- 



H 



eff 



■A< 



JV- 



"l)/2 V + 
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(36) 

Here we have not included the renormalization fac- 
tor dYi/duj in eq. (36), because at T = it does 
not affect the dc conductance and charge displacement 
N c \. The eigenvalue e* of H^ satisfies the condition 
det {G qp (e*)} _1 = 0. Owing to the inversion symmetry, 
this equation can be factorized, as eq. (21) 

H [det IT {u N ) - v 2 k( N -VI 2 g N {e*) det ^(wjv) 



0, 



(37) 



where the product runs over the two partial waves, 7 = 
"even" and "odd". In eq. (37), the argument for the de- 
terminants is ujn = £* A - ^ -1 )/ 2 , which becomes zero in 
the limit of N — > 00. The Green's function g N (cu) is de- 
fined in the site n — at the interface of an isolated lead 
consisting of N + 1 sites, and it can be expressed in the 
form g N {u) = Em=o \r°m(0)\ 2 / (w - e m ), where e m and 
(p m {n) are the eigenvalue and eigenfunction for the iso- 
lated lead. Furthermore, using eq. (37), the eigenstates 
of H™ can be classified according to the parity. Thus, 
the two parameters K C von and K dd, which are defined in 
eq. (22), can be deduced from the low- lying eigenvalues 
for the quasi-particles of the even and odd parities, £o Von 
and e* dd , as 



log A, (33) ^ 



det K 7 (0) 
Tdet Kjjo) 



= -__ lim D\( N ~ 

YD N^oo 



1)11 



(38) 



for 7 = "even" and "odd". Note that v 2 /(TD) = 2A A /n 
from the definition eq. (33). With these eigenstates, the 
quasi-particlc Hamiltonian is written in a diagonal form 



H (N) 

qp 



EE 



'-h,lj 



9 ' 



(39) 

where E g is the ground-state energy, and e*j (£/j.; 7 ) 
is the l-th excitation energy of a single-particle (single- 



J. Phys. Soc. Jpn. 



Full Paper 



A. Ogupj, Y. Nisikawa, and A. C. Hewson 



5 



V t t V 



Fig. 4. Schematic picture of a series connection for Nq 



hole) state for parity 7. Thus, the parameters K C von and 
Kodd can be calculated substituting the NRG result of a 
low-energy eigenvalue, e* , or — e* h , , into the right-hand 



side of eq. (38). Then, the conductance, local charge, and 
phase shifts <5 CV cn and <5 dd, can be calculated using the 
equations given in §2. 

Alternatively, one can calculate the conductance di- 
rectly from the current-current correlation function with- 
out using the fixed-point Hamiltonian of the local Fermi 
liquid. 11,12 It is more general, and is applicable also at 
finite temperatures. For the systems consisting of a num- 
ber of interacting sites Nc, however, the number of eigen- 
states needed for carrying out the NRG iterations be- 
comes large, and the calculations of the expectation val- 
ues require a rather long computer-time. Therefore, for 
studying the ground-state properties, the formulation de- 
scribed above is much easier. Moreover, our formulation 
requires only the eigenvalues to deduce the phase shifts 
via eq. (38). This is more efficient than deducing them 
from the expectation values of the local charge {rid^o), 
which would have to be determined by an additional cal- 
culation. 

4. Phase shifts for triple-dot systems 

We now apply the formulation described in the above 
to a finite Hubbard chain coupled to two noninteracting 
leads. Specifically, we concentrate on a triple dot as il- 
lustrated in Fig. 4. The explicit form of the interaction 



Hamiltonian is Hq 



d,i-\ n d,il 



with Nc 



We assume that the bare hopping matrix elements ty, 
which are defined in eq. (2) for Tt c , are described by the 
nearest-neighbor hopping t, and assume that all other 
off-diagonal elements to be zero. The onsite energy of the 
interacting cites is given by — tu = td for 1 < i < Nc, 
and the origin of the energy is chosen to be \i = 0. 

4-.1 Series connection 

We first of all discuss the noninteracting case U = 0. 
In Fig. 5, the conductance Series, local charge in the 
triple dot N e \, and the phase shifts are plotted as func- 
tions of td/t for (solid line) T/t = 0.12 and (dashed line) 
r/t = 0.5. The conductance shows sharp peaks when the 
resonance states in the triple dot cross the Fermi level. 
Among the three conductance peaks, the one in the mid- 
dle corresponds to an orbitals with the odd parity, and 
the remaining two belong to the orbitals with the even 
parity. For U — 0, the two phase shifts are given by 
cot Seven = (e 2 ; _ 2f2 ) / ( Te d), and cot S odd = e d /T. The 
phase shifts show a step of the height 71", when the reso- 
nance state of the corresponding parity crosses the Fermi 
level. Simultaneously, the local charge N e \ also shows a 
staircase behavior, and the step corresponding to two 
electrons, spin up and spin down, occupying the reso- 
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Fig. 5. The conductance g se ries, and total charge in the triple- 
dot N c i in the noninteracting case U = as functions of e^/t for 
(solid line) T/t = 0.12 and (dashed line) T/t = 0.5. The phase 
shifts 26 cvcii/tt and 2S odd /n are also shown in (b). 



nance states. For large T, the resonance peaks become 
broad, and the steps in the phase shifts vanish gradually. 

For interacting case U 7^ 0, we have carried out the 
NRG iterations by introducing the bonding and anti- 
bonding orbitals also for the leads; (/^ Ra ± p n Lcr )/V2- 
In the present study, instead of adding these two orbitals 
at n = N + 1 simultaneously for constructing the Hilbert 
space for the next NRG step, we add the bonding or- 
bitals first and retain typically 3600 low-energy states 
after carrying out the diagonalization. Then, we add the 
remaining anti-boding orbitals, and carry out the trun- 
cation again keeping the lowest 3600 eigenstates. This 
truncation procedure preserves the inversion symmetry, 
and in the case of odd Nc it also preserves the particle- 
hole symmetry at = —U/2. With this procedure, the 
discretized Hamiltonian Hn can be diagonalized exactly 
up to N — 1, where the total number of the sites in the 
cluster is N c + 2(N + 1) = 7 for the triple dot N c = 3 
as illustrated in Fig. 3. For N > 2, we have confirmed 
that the eigenvalues in the noninteracting limit U — > 
are reproduced sufficiently well with this procedure for 
the discretization parameter of the value A = 6.0. 

We have also confirmed that the fixed-point eigenval- 
ues of the NRG for the triple dot can be described by 
the quasi-particle Hamiltonian defined in eq. (39) for all 
parameter values we have examined. It justifies the as- 
sumption of the local Fermi liquid we have made to derive 
eq. (38). In Fig. 6, the NRG results for (a) the conduc- 
tance ff ser j es , (b) phase shifts 5 ovon , 5 dd and total charge 
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Fig. 6. The conductance g series , local charge N c \, phase shifts 
even/it and 2S ^/tt for the triple dot as a function of €d/U 
for U/(2nt) = 1.0 and T/t = 0.12. The dashed vertical lines in 
(b) show the values of at which N e \ jumps in the "molecule 
limit" T = 0. For NRG, we use t/D = 0.1 and A = 6.0. 
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Fig. 7. The conductance 9 aeries , local charge N e i, phase shifts 
2S evcu /n and 2(5 dd/7r for the triple dot as a function of td/U 
for U/(2nt) = 1.0 and T/t = 0.5. The dashed vertical lines in (b) 
show the values of at which 7V i jumps in the limit of T = 0. 
For NRG, we use t/D = 0.02 and A = 6.0. 



in the triple dot N e \ are plotted as a function of e<j/Z7 
for U/(2nt) = 1.0. Here the hybridization is relatively 
small T/t — 0.12 compared to the hopping matrix el- 
ement t, which is taken to be t/D = 0.1 in the NRG 
iterations. The vertical dotted lines in (b) correspond to 
the values of td at which N e \ jumps discontinuously in 
a 'molecule' limit r — > 0. The local charge N c \ shows 
a step behavior near the dotted lines. Specifically, due 
to the Coulomb interaction, the steps emerge also for 
odd N c \. These odd steps reflect the tt/2 steps of the 
phase shifts. The conductance shows the typical Kondo 
plateaus of the Unitary limit g ~ 2e 2 /h in the regions 
of €d corresponding to the odd occupancies N c \ ~ 1, 3, 
5. In contrast, the conductance shows wide minima for 
even occupancies N c \ ~ 2, 4. These features are linked to 
the behavior of the phase shifts; <5 C vcn — ^odd — tt/2 for 
odd N c \, and 5 ovcn — (5 oc jd — or ir for even iV e i. 16 

Among the three conductance plateaus, the one in the 
middle is wider than the others. This seems to be caused 
by the suppression of the charge fluctuation near half- 
filling. It is expected that the plateau becomes wider for 
some particular values of N e \, at which the charging en- 
ergy defined with respect to the unconnected limit T = 
is large. Thus, a wide plateau may emerge for the fill- 
ing near the metal-insulator transition, which generally 
depends on the spatial range of the repulsive interaction. 

We next consider the case where the triple dot is cou- 



pled strongly to the leads via a large hybridization. In 
Fig. 7, the NRG results obtained at T/t = 0.5 are shown, 
where the value of the Coulomb interaction is unchanged 
U / (2irt) = 1.0. The local charge JV e i, which shows a stair- 
case behavior for small T in Fig. 6, becomes now a gentle 
slope in Fig. 7 (b) due to the large hybridization. Cor- 
respondingly, the shoulders of the conductance plateaus 
become round, and the valleys become shallow. Never- 
theless, near half-filling, the phase shift for the odd par- 
tial wave still shows a weak tt/2 step, i.e., <5 dd — 7r /2 near 
half-filling. It keeps the conductance peak at the center 
flat. Note that in this calculation the hopping matrix el- 
ement is taken to be t/D — 0.02, but the change in this 
ratio does not affects the low-energy properties so much. 

We also examine the case with a weak interaction 
U/{2nt) = 0.2 and small hybridization T/t = 0.12. The 
results are shown in Fig. 8. It shows an early stage of a 
development of the Kondo plateaus. In the upper panel 
(a), the dashed line corresponds to the conductance in 
the noninteracting case: for comparison it is shifted to 
the negative energy region by — U/2 so that the positions 
of the central peak coincides with that for the interact- 
ing case, and also the energy is scaled by U using the 
value U/(2itt) = 0.2. Due to the Coulomb interaction, 
the conductance peaks become wider than that for the 
noninteracting electrons, and the form of the peaks de- 
viates from the simple Lorentzian shape. Furthermore, 
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Fig. 8. The conductance <7 8eries , local charge N c \, phase shifts 
2<5even/i" and 2<5 0< j < j/"" for the triple dot as a function of e^/U 
for U/(2nt) = 0.2 and Y/t = 0.12. In (a) the dashed line cor- 
responds to the results for U = 0: it is shifted to the negative 
energy region by —U/2 so that the position of the central peak 
coincides with that of the interacting case, and also the energy 
is scaled by U using U/(2nt) = 0.2. The dashed vertical lines 
in (b) show the values of at which jV e i jumps in the limit of 
T = 0. For NRG, we use t/D = 0.1 and A = 6.0. 



the separation of the peaks increases. Since U is not so 
large in this case, the local charge N e \ becomes flat only 
for even occupancies N e \ = 0,2,4,6. It reflects the be- 
havior of the phase shifts, i.e., both 5 cvcn and <5 dd do 
not show the clear ir/2 steps. Nevertheless, the slopes in 
between the tt steps become longer compared to the ones 
in the noninteracting case, and weak precursors for the 
development of the ir/2 steps can be seen in these slopes. 

4-2 Parallel conductance 

The parallel conductance through the tripled dot as 
illustrated in Fig. 9, can be evaluated via eq. (29) with 
the phase shifts for the series connection given in the 
above. This is because, as mentioned in §2, the matrix 
Green's function G(z) for the interacting sites for the 
parallel connection coincides with that for the series con- 
nection owing to the inversion symmetry for the leads in 
the horizontal direction, as described in Fig. 2. Thus a 
site-diagonal quantity such as local charge N e \ becomes 
the same with that for the series connection as mentioned 
in §2. The conductances, however, are different. It opens 
a way to determine the values the phase shifts <5 e ven and 
<5odd from measurements of the series and parallel con- 



V/y/2 V/V2 



V/V2^V/V2 



Fig. 9. Schematic picture of a parallel connection for No = 3. 



ductances using eqs. (27) and (29). 

We first of all consider the noninteracting case. In Fig. 
10, the parallel (solid line) and series (dashed line) con- 
ductances for L7 = are plotted as functions of e/t for 

(a) Y/t = 0.12 and (b) Y/t = 0.5. In (a), the two curves 
almost overlap each other, and the difference between the 
parallel and series conductances can be seen only near the 
two valleys corresponding to even JV e i. The both conduc- 
tances show a similar feature due to the three resonance 
states. At half-filling, Cd = for noninteracting electrons, 
the two curves contact with each other. This is because 
the phase shifts take the values <5 OV cn = ^ and 5 dd = tt/2. 
The phase shifts do not show such an exact synchronism 
at the side resonance peaks corresponding to N c \ ~ 1 and 
5, and there the parallel conductance become larger than 
2e 2 /h, which is possible because two conducting channels 
contribute to the current in the parallel connection. The 
parallel conductance becomes larger than the series con- 
ductance for \ed\ > 2t. This is because the both phase 
shifts show the same asymptotic behavior, 5 cvcn ~ Y/ej, 
and J dd — r/ed, for — > oo, and a similar cancellation 
occurs for the series conductance also in the opposite 
limit td — > — oo. The difference between the parallel and 
series conductances becomes small for small T. 

In Fig. 11, the (solid circle) parallel and (open square) 
series conductances in the interacting case U /{2nt) = 1.0 
are plotted as a function of td/U for (a) Y/t = 0.12 and 

(b) Y/t = 0.5. The difference between g pa raiici and emeries 
for interacting electrons is similar, qualitatively, to that 
in the noninteracting case, although the appearance of 
the Kondo plateaus changes the overall feature of the Cd 
dependence. In (a) the difference are visible only near the 
conductance valleys. For large Y, however, the difference 
becomes larger, and it can be seen whole region plotted 
in (b) except near the central plateau. As in the nonin- 
teracting case, the parallel conductance becomes larger 
than 2e 2 /h at the side Kondo plateaus corresponding to 
N c \ ~ 1 and 3 due to the contributions of the two con- 
ducting modes. The conductance valleys corresponding 
to even N e \ become deeper for the parallel conductance 
than that of the series conductance. However, the differ- 
ence is small for small Y. 

5. Discussion 

The results presented in the above have been obtained 
at zero temperature. Nevertheless the Kondo behavior we 
have discussed can be seen at low enough temperatures 
T < Tk, although the height of the Kondo plateaus will 
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Fig. 10. The conductance in the parallel (solid line) and series 
(dashed line) connections for noninteracting electrons U = 
as functions of e^/t, where (a) Y/t = 0.12 and (b) Y/t = 0.5. 
Note that two conducting channels contribute to the parallel 
conductance, so that the upper bound for 5p ara iiel is 

not 2e 2 /h. 
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Fig. 11. The conductance in the parallel (solid line) and series 
(dashed line) connections for interacting electrons U/(2nt) = 1.0 
as functions of e d /U, where (a) Y/t = 0.12 and (b) Y/t = 0.5. 
The two arrows in (a) indicate the points, for which the energy 
flow is described in Fig. 12. 



decrease gradually with increasing T from the Unitary- 
limit value as that in the single impurity case. 12 At 
T ^> Tk, the plateaus finally disappear and the con- 
ductance shows the oscillatory behavior of the Coulomb 
blockade with the peaks, the height of which is about 
half of the Unitary- limit value. 16 The Kondo tempera- 
ture Tk decreases with increasing U, and a large num- 
ber of the NRG iteration N is required to get to the 
Fermi-liquid fixed point that determines the low-energy 
Kondo behavior for large U. 13 The value of Tk should be 
different for different plateaus. Specifically, for the cen- 
tral plateaus near half-filling, Tk becomes small for the 
Hubbard chain with large size Nq because the develop- 
ing Mott-Hubbard (pseudo) gap disturbs the screening 
of the local moment. 

In order to estimate the Kondo energy scale for dif- 
ferent plateaus, the flow of the low-lying eigenvalues of 
Hn/D for the triple-dot is plotted in Fig. 12 as a func- 
tion of odd N for U/(2irt) = 1.0 and T/t = 0.12. This 
parameter set is the same as the one used for Fig. 6 and 
Fig. 11 (a). The figure 12 (a) shows the flow at the central 
Kondo plateau for ed/U = —0.5, and the lower panel (b) 
shows the flow at the right plateau for ed/U ~ 0.15625. 
The positions corresponding to these two values of Cd are 
indicated in Fig. 11 (a) by the arrows. In the figure 12, 
the crossover from the high-energy region to low-energy 
Fermi-liquid regime can be seen clearly. For large N, the 



energy levels converge to the fixed-point values, and the 
many-body eigenvalues of Hn / D have one-to-one corre- 
spondence with the free quasi-particle excitations. 9 ' 10 ' 13 
The number of NRG iterations N* needed to reach the 
low-energy regime in each of the plateaus is estimated 
to be (a) N* ~ 30, (b) N* ~ 15, and N* ~ 15 also 
for the left plateau, for which the energy flow is identi- 
cal to Fig. 12 (b) except for the values of some quan- 
tum numbers assigned to the lines. The crossover occurs 
at an energy T* ~ DA~( N _1 )/ 2 , where a factor of the 
form A^^ 1 )/ 2 emerges to recover the energy scale of 
the original Hamiltonian TL from the discretized version 
Hm defined in eq. (30). The crossover energy reflects the 
width of the Kondo resonance, so that T* ~ Tk apart 
form a numerical factor of order 0(1) which depends 
on the precise definition of Tk- Consequently, these re- 
sults show that the Kondo energy scale for the central 
plateau at half-filling is much smaller than that for the 
side plateaus away from half-filling. 

In order to calculate the conductance for whole tem- 
perature ranges, one generally has to use several different 
approaches. The behavior at high temperatures T > Tk 
can be captured, for instance, by the exact diagonal- 
ization of small clusters, such as the one examined by 
several authors. 24 ' 25 At low temperatures, however, one 
needs the information about the low-energy excitations. 
From the Fermi-liquid regime to the crossover region near 
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Fig. 12. Low energy levels of Hjy/D are plotted as a function of 
odd N for (a) e d /U = -0.5 where N el = 3, and (b) e d /U = 
0.15625 where iV e i ~ 1. Other parameters are the same with 
those in the case of Fig. 6: U/{2irt) = 1.0 T/t = 0.12, t/D = 0.1 
and A = 6.0. The label (Q, 2S, p) corresponds to the charge Q, 
spin S, and parity p (= even, or odd). 



Tk, the conductance can be calculated from the current- 
current correlation function using NRG method. 11 ' 12 
The ground-state value of the conductance can be cal- 
culated from the phase shifts, which can be deduced ef- 
ficiently from the low-energy eigenvalues using eq. (38). 

In summary, we have described the method to deduce 
the phase shifts for finite interacting systems attached to 
nonintcracting leads from the fixed point Hamiltonian of 
NRG. Our approach assumes the inversion symmetry, 
but does not assume a specific form of the Hamiltonian 
nor electron-hole symmetry, and thus it is applicable to a 
wide class of quantum impurities. We apply the method 
to a triple quantum dot of Nc = 3 connected to two non- 
interacting leads, and calculate the dc conductance away 
from half-filling as a function of the onsite energy e,j. At 



T = 0, the conductance shows the Kondo plateaus of the 
Unitary limit g ~ 2e 2 /h at the values of gate voltage 
corresponding to odd numbers of electrons AT i, while it 
shows the wide minima for even N c \. It seems to be natu- 
ral to expect that these low-temperature properties seen 
at T < Tk are common to the Hubbard chain of finite 
size Nc attached to Fermi-liquid reservoirs. 
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